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The dynamic susceptibility of propylene carbonate in the moderately viscous regime above Tc 
is reinvestigated by incoherent neutron and depolarised light scattering, and compared to dielectric 
loss and solvation response. Depending on the strength of a relaxation, a more or less extended (3 
scaling regime is found. Mode-coupling fits yield consistently A — 0.72 and Tc = 182 K, although 
different positions of the susceptibility minimum indicate that not all observables have reached the 
universal asymptotics. 



I. INTRODUCTION 

A. Motivation 

The glass transition, an essentially dynamic phe- 
nomenon, can be described as the slowing down and even- 
tual freezing of a relaxation. According to mode-coupling 
(MC) theory the physical origin of this process must 
be sought on a picosecond scale where long-ranged trans- 
port starts to evolve from vibrational short-time dynam- 
ics. 

In the long-time limit, MC theory reproduces the well- 
established phenomenology of a relaxation. New results 
are obtained for shorter times. In particular, the theory 
predicts a change of transport mechanism around a cross- 
over temperature Tc, located in the moderately viscous 
liquid phase well above the conventional glass transition 
temperature Tg. On cooling towards Tc, particles spend 
more and more time being trapped in transient cages; 
this process, labelled fast P relaxation, is predicted to 
obey remarkably universal scaling. 

In a couple of structural glass formers, MC predic- 
tions have been confirmed |^,^ primarily by different 
scattering techniques More recently, the GHz- 

THz dynamics also became accessible by dielectric spec- 
troscopy 0,^. Although results are in accord with the 
MC scenario, in several cases the data did not fall to- 
gether with dynamic susceptibilities from scattering ex- 
periments [^,0 which is in conflict with the asymptotic 
factorisation property of MC /3 relaxation. Recent theo- 
retical developments suggest explanations on the basis of 
corrections to scaling ||ll|,|l2| and orientational degrees of 
freedom For experimental tests, more detailed 

comparisons of different observables are needed. 

In this context, we performed incoherent neutron and 



depolarised light scattering experiments on propylene 
carbonate (PC, 4-methyl-l,3-dioxolan-2-on, C4O2H6), a 
fragile glass former (Tg — 160 K) with low molecular 
weight {M — 102.1) which has already been studied by 
various experimental techniques. A synopsis of available 
data shall be given in subsection I.C after collecting the 
essential results of MC theory in LB. 

B. The Mode-Coupling Cross-Over 

In its simplest ("idealised") formulation, MC theory 
describes an ergodic-to-nonergodic transition at Tc. On 
the low-temperature side, the onset of fast /? relaxation 
leads to an anomalous decrease of the Debye- Waller or 
Lamb-Mofibauer factor, 

f, = f^ + K\a\'/\a>0, (1) 

with a reduced temperature a — (Tc — T)/Tc. On the 
high-temperature side, the time constant of a relaxation 
diverges with a fractal exponent 7, 

r cx IctI"'' , fj < 0. (2) 

Such a transition has actually been observed in a col- 
loidal suspension |]l6| . In a structural glass former, the 
singularities (|^) and (||) are smeared out because acti- 
vated hopping processes restore ergodicity ||l^. Under 
these limitations, integral quantity like fq or Tq do not 
allow for decisive tests of theory. One rather has to study 
the full dynamics, as represented by the dynamic suscep- 
tibility x^'(cij), or any other dynamic variable coupling to 
it. 

Stronger predictions are made for the fast (3 regime: 
around the minimum between a peak and vibrational 
excitations, in a temperature range close enough to, but 
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sufRciently above Tc, any susceptibility is expected to 
reach the same asymptotic hmit 



II. NEUTRON SCATTERING 
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A. Experiments 



where the scahng function gx is fuUy determined by one 
single parameter A fl^ ]. Amplitude and frequency scale 
of (y) should become singular on cooling towards Tc, 

X<T oc |(T|l/^ and w,, cx |(T|i/2a^ ^ ^ Q 

where the exponent a, just as 7 in (|^), is determined 
by A. 

C. Previous Studies of Propylene Carbonate 

Table 1 summarises previous studies of fast dynamics 
in PC. All authors reported at least partial accord with 
MC predictions. However, not all MC interpretations 
were consistent with each other. It is an essential result 
of idealised MC theory that in the asymptotic (3 regime 
all dynamic observables show the same spectral distribu- 
tion and the same temperature dependence. Therefore, 
one material is characterised by just one A and one T^. 
Some of the Tc reported for PC are therefore incompati- 
ble with MC theory. 

As with other materials, early scattering experiments 
[ p9| , p0t concentrated on the square-root singularity of fq. 
As in other materials, this singularity remains elusive: in 
Brillouin scattering, data fitting depends on uncontrolled 
approximations for the memory function Similarly, 
in neutron scattering a determination of Tc from (Q) 
works at best if the full lineshape on the cr > side is 
known from high-resolution spectroscopy p2[ . Therefore, 
the isolated results Tc = 210 or 270 K can be discarded 
from further consideration. 

In comparison, the determination of Tc from viscosity 
or a-relaxation data works better. However, available 
data do not allow for an independent determination of A, 
and even for one given value A = 0.70, results vary be- 
tween Tc = 180 and 196 K. 

A dynamic susceptibility has been measured first by 
depolarised light scattering yielding A = 0.78 and 
Tc = 187. These values have shown to be consistent with 
dielectric loss spectroscopy [p|, p3|j2^ ]. However, a time- 
domain optical measurement of the solvation response of 
a solute molecule found a significantly lower Tc if the 
light scattering value of A = 0.78 was assumed p5| , p6t . 

In the following section II, we study the fast relaxation 
regime by neutron scattering. The A and Tc we obtain 
differ significantly from results of the other dynamic mea- 
surements which motivates us to remeasure some light 
scattering spectra (sect. HI) and to reanalyse dielectric- 
loss and solvation- response data (sect. IV). Only then, 
all susceptibilities will be compared in section V. 



Inelastic, incoherent neutron scattering has been mea- 
sured on the time-of-flight spectrometers Mibemol at the 
Laboratoire Leon Brillouin, Saclay, and Neat at the 
Hahn-Meitner-Institut, Berlin. 

On Mibemol, the counter-rotating choppers were op- 
erated with 10000 rpm. With an incident neutron wave- 
length of Ai = 8.5 A, we obtained a resolution (fwhm) of 
6.5 — 9.5 GHz, depending on angle. At the Berlin reac- 
tor, the flux delivered by the undermoderated cold source 
decreases rather fast for wavelengths beyond 5 or 6 A; 
therefore Ai = 5.5 A was chosen. With counter-rotating 
choppers at 20 000 rpm, we achieved nevertheless a res- 
olution of 11 — 13.5 GHz. Total count rates were of the 
same order for both instruments; a precise comparison 
cannot be made because different sample geometries were 
used. 

For the Mibemol experiment, we used an AI hollow 
cylinder |p7|| with 30 mm outer diameter and a sam- 
ple layer of about 0.1mm thickness. In this container, 
the sample crystallised partially during a 190 K scan; at 
175 K, as far as one can tell from the elastic structure 
factor of a predominantly incoherent scatterer, complete 
or nearly complete crystallisation occured within 2 hours. 
On the other hand, after rapidly cooling from 260 K, it 
was possible to perform a 2 K scan in the fully amorphous 
state. 

On Neat, we tried to remeasure the dynamics of the 
supercooled state below 200 K. To this purpose, we filled 
the sample into nearly 200 thin capillaries (soda lime, in- 
ner diameter 0.2 mm, Hilgenberg). As in other liquids, 
this packaging proved highly successful; comparison to a 
vanadium scattering and visual inspection showed that 
down to 168 K no crystallisation occured. Because the 
detectors arc located at only 250 cm from the sample, the 
resolution of Neat is very sensitive to flightpath differ- 
ences. We therefore abstained from using a hollow cylin- 
drical geometry; instead, we placed the capillaries in a 
rectangular holder with a 30 mm base length which was 
then mounted at 45° with respect to the incoming beam. 
For both experiments, the sample material, 1,2-propylene 
carbonate (99.7%, Sigma- Aldrich) , was loaded under an 
inert gas. 

After converting the raw data to S{29,uj), they were 
binned in about 30 angular groups; a non-equidistant fre- 
quency binning was imposed by requiring statistic fluctu- 
ations to fall below a given mark. Only then was the con- 
tainer scattering subtracted. The absolute intensity scale 
was taken from the elastic scattering of the Mibemol 2 K 
scan; for Neat, detectors were first calibrated to vana- 
dium before the overall scale was fit to Mibemol at 210 K. 
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The scattering law S{29,llj) depends still on Aj. Only 
after interpolation to constant wavenumbers q one ob- 
tains spectra S{q,uj) which are independent of the kine- 
matics of scattering, and only then a direct compar- 
ison between our two experiments becomes possible. 
This comparison is performed explicitely at 207/210 K in 
Fig. |l|. In the quasielastic scattering range, up to some 
100 GHz, the accord between Mibemol and Neat is ex- 
cellent. 

B. The Factorisation Property 

For analysing the data, especially in the quasielastic 
range, it is advantageous to visualise the scattering law 
as a susceptibility 

Xqiuj) ^ S{q,u;)/n{Lj) (5) 

with the Bose factor n{uj) = {exp{hLLj / kBT) — 1) ^. 

Around the susceptibility minimum, mode coupling 
theory predicts that the (3 line shape is asymptotically 
the same for all experimental observables which couple 
to density fluctuations; thus a fortiori it must be the 
same for incoherent neutron susceptibilities measured at 
different wavenumbers q. We therefore expect Xqi^) 
to factorise into an w-independent amplitude and a q- 
independent spectral function, 

= (6) 

On different theoretical grounds, such a factorisation is 
also expected for incoherent scattering from harmonic vi- 
brations in lowest order of mass expansion |2^ , |2^ . 

There are several ways to test (^) and to determine hq: 
e. g., one may iteratively construct a model for and 
fit it to the individual Xq{^) data sets. More simply, hq 
can be calculated from a least-squares match between 
neighbouring q cuts [|9|. For the present data, we find 
that hq does not depend on the chosen procedure, nor 
does it vary with the frequency subrange from which it 
is determined. 

A first surprising result is then the strictly linear 
wavenumber dependence of hq (Fig. |^). Within MC 
theory, q dependences can be calculated only if a specific 
microscopic structure is put in. Within harmonic the- 
ory, however, there is a clear expectation that hq (x q^. 
Astonishingly enough, this dependence is never seen. 
Instead, in at least two other molecular systems a lin- 
ear hq (X q dependence is found though the temperature 
dependence of THz modes indicates pure harmonic be- 
haviour ||^,|l|- From other time-of-fiight studies [||J32) 
we suspect and a simulation confirms that this is 
mainly a multiple-scattering effect. Anyway, the linear- 
ity is so accurate that in our analysis we shall use (|^) with 
hq = q/qo (with an arbitrary normalisation qo — 1 A^^) 
instead of employing the empirical values. 



Fig. ^ shows that the factorisation holds over the full 
experimental data range, at least from 30 to 2500 GHz, 
except for the onset of a relaxation at the highest tem- 
perature (which is associated with long-ranged transport 
and depends therefore strongly on q). This allows us 
to condense our huge two-dimensional data sets into q- 
independent functions x"i^) with much improved statis- 
tics; only on this basis is it possible to proceed with a 
quantitative analysis of noisy data accumulated in rela- 
tively short scans. 

Before proceeding with the analysis we have to note 
that the x"('^) determined by this method and shown 
for 207/210 K in Fig. ||b reveals considerable discrepan- 
cies between Mibemol and Neat at frequencies above 
about 400 GHz. Although these deviations passed unno- 
ticed on the double logarithmic scale of Fig. |l|b, they are 
already present in the raw data. While we can exclude 
container scattering and dark counts as possible causes, 
two other frequency-dependent effects are likely to con- 
tribute: multiple scattering [^3| and an inaccurate detec- 
tor efficiency correction [Q. Multiple scattering cannot 
be corrected for unless a comprehensive theoretical model 
of S{q,uj) is used as input; therefore, it presents a funda- 
mental limitation to the determination of spectral shapes 
by quasielastic neutron scattering. 



C. Master Curves for f3 Relaxation 

The g-independent x"i^) can now be used to test the 
scaling form (||). We use fits with a fixed A to deter- 
mine Xcr and Wct- With these values, a master curve 
x!' (yi I Lo I Xa is constructed, and from the master curve 
the scaling range is read off which is then used for im- 
proved fits to the original data. Figure ^ shows fits to the 
original x"('^) for A — 0.72, and Figure ^ shows master 
curves for three A's. 

In principle, a self-consistent A can be determined in 
an iterative procedure from free fits to the master curves. 
For our 220 — 260 K data, A tends towards values around 
0.69 (Figure ||a), but the convergence is erratic, and the 
outcome depends on the subjective decision of which 
points to include in the master curve. 

With the data on hand, a more restrictive determina- 
tion of A is possible from the temperature dependence (Q). 
The insets in Figure ^ show that a consistent linear be- 
haviour Xa^ cx; Wo-^" oc (7 is found only for A = 0.72 and 
Tc = 182 K; the exponent a — 0.32, which corresponds to 
A = 0.72, is confirmed by cross-checking \T^Xa vs. Iulu^ 
(inset of Figure |4|) . 

Figure |^ demonstrates that the value A = 0.78 sug- 
gested by preceeding light scattering and dielectric loss 
measurements do not give a good description of the neu- 
tron scattering data: the master curve is of poorer qual- 
ity than for A = 0.69 or 0.72, especially at frequencies 
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around and below the minimum; free fits with g\ show 
a clear trend towards smaller values of A; and the scales 
Xa and Wo- do not consistently follow (||). 

This discrepancy motivates us to remeasure some light 
scattering spectra around the susceptibility minimum. 

III. LIGHT SCATTERING 
A. Previous Studies 

Light scattering, just as neutron scattering, comprises 
coherent and incoherent contributions. Since the wave- 
length of visible light is much longer than molecular 
dimensions, coherent scattering arises only from sound 
modes, giving rise to discrete Brillouin lines. Incoher- 
ent scattering, on the other hand, sees local motion and 
yields a continuous spectrum. In a very first approxima- 
tion, this spectrum can be interpreted as if it were a q 
averaged neutron scattering law; microscopic models sug- 
gest that the scattering mechanism involves four-point 
density correlations and/or rotational motion. 

Brillouin scattering yields the velocity and damping 
constant of sound waves. In principle |]35[|, the sound dis- 
persion through the glass transition reveals the strength 
of the a relaxation so that it is possible to read off the 
Debye- Waller factor fq from fits to the Brillouin lines. 
In practice, this led to an inconsistent estimate of Tc in 
PC [^; for a reliable determination of fq{T), one should 
not only determine one limiting sound velocity by ul- 
trasonic experiments, but also provide an independently 
determined memory function as input pl[ |. 

More direct information on the microscopic dynamics 
is obtained from the incoherent continuous spectra. In 
order to suppress the much stronger Brillouin lines, these 
spectra are preferentially gathered in depolarised (VH) 
back-scattering geometry. In an extensive study of PC, 
besides a thorough discussion of Brillouin scattering, Du 
et al. have also measured VH spectra over a wide 
frequency band, made accessible by combining a tandem 
interferometer with a grating monochromator. The mea- 
surements extended over a large temperature range, in- 
cluding a relaxation as well as hopping processes below 
Tc, and were analysed within both the "idealised" and 
"extended" version of MC theory. 

While these measurements gave a broad overview of 
the dynamics of glass-forming PC, they did not yield 
spectral lineshapes with the precision we need now for a 
quantitative comparison with neutron scattering results. 
In particular, the exponent parameter A — 0.78, obtained 
within idealised theory from a global fit to T > Tc data, 
was given with a relatively large error range of ±0.05. 
Furthermore, as in other broad band light scattering 
studies performed until quite recently, the tandem in- 
terferometer was used in series with an insufficient band- 
pass which did not fully suppress higher-order transmis- 



sions of the interferometer | |3q , |37| . For these reasons, we 
remeasured some depolarised spectra, concentrating on 
the temperature and frequency range of the asymptotic 
(3 regime. 

B. Experimental Set-up 

Experiments were performed in Garching on a Fabry- 
Perot six-pass tandem interferometer. The instrument, 
bought from J. R. Sandercock, has been modified in sev- 
eral details to allow for stable operation and high con- 
trast. The six-pass optics has been placed in a thermally 
isolating housing, and the scanning stage is actively tem- 
perature stabilised. The analogue-electronic stabilisation 
of the interferometer piezos is replaced by computer con- 
trol. Entrance and exit pinholes are spatially separated 
from the six-pass optics, so that the most critical aligne- 
ments can be done without disturbing the interferometer 
operation. By placing additional masks in the six-pass 
optics, in particular on mirror surfaces, the cross-talk 
between different passes could be reduced by several or- 
ders, and a straylight rejection of better than 10^^ was 
achieved. 

Depending on the free spectral range, the instrument is 
used in series with an interference filter of 0.15 or 1 GHz 
band width (Andover); these filters are placed in a spe- 
cial housing with active temperature stabilisation. Fur- 
thermore, to account for long-term drift, the instrument 
function is redetermined periodically by automatic white- 
light scans. 

Although the spectrometer guaranteed excellent stray- 
light rejection, special care was taken to prevent direct or 
diffuse reflections of laser light from entering the instru- 
ment. Therefore, instead of 180° back-scattering, a 169° 
VH geometry was chosen. From the intensity transmit- 
ted through the "ghosts" of the tandem instrument, we 
conclude that the straylight was about 10^ times weaker 
than the inelastic scattering from PC, and completely 
negligible compared to the detector's dark count rate of 
about 2.5 sec^^. After subtraction of these dark counts, 
and normalisation to the corresponding white light scans, 
the resulting spectra showed excellent detailed balance 
symmetry. 

The sample material was from the same source as 
for the neutron scattering experiment, and was vacuum 
sealed in a Duran cell. 

C. Susceptibilities Around ujcr 

For a precise determination of the spectral lineshape, 
subsequent measurements over different spectral ranges 
were performed after stabilising the temperature over 
night. The most restrictive determination of the expo- 
nent parameter was possible at T = 216 K. As shown 
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in Fig. |[ the matching of the three overlapping spec- 
tral ranges is excellent. Except for the leaking VV Bril- 
louin mode, the VH susceptibility is described over more 
than two decades by the mode coupling asymptote ^a- 
Fits yield A = 0.72 ± 0.01, which is confirmed by mea- 
surements at other temperatures as well as by a bulk of 
earlier experiments we had performed under less ideal 
experimental conditions. This result is at the margin of 
the error range in the literature value 0.78 ±0.05 the 
figure shows that 0.78 itself is clearly incompatible with 
our present data. 

Secondly, for two different spectral ranges we measured 
temperature series. Fig. |^ shows some of the composite 
susceptibilities. From fits with fixed A = 0.72, the fre- 
quencies oja were obtained. Then, the individual suscep- 
tibilities, measured around interferometer mirror spac- 
ings zq — 0.8 and 2.4 mm were fitted with fixed A and lo^ 
so that we obtained two independent data sets for the 
amplitudes Xo-- 

The temperature dependence of Wo- is shown in Fig. |^a. 
Data between 190 and 230 K extrapolate to the same 
Tc = 182 ± 1 K as found from neutron scattering. Ex- 
tending the fit range to higher temperatures leads to a 
Tc which is about 2K lower. Again, these results are 
marginally compatible with Tc — 187 ±5K from the ear- 
lier light scattering study pl| ]. 

Measuring amplitudes in light scattering is difficult, 
and the prediction (|^) can be verified only over a reduced 
temperature range. Results are visualised best in a log- 
arithmic plot of Xct vs. T — Tc (Fig. ||b) which suggests 
there are two different regions in which Xir is proportional 
to |crp/^, separated by some step. Without further ex- 
periments, we must leave open whether this step comes 
from the sample or presents an experimental artifact, due 
for instance to distortions of the optical paths. We note, 
however, that after 30 hours and a full temperature cycle, 
we were able to reproduce an amplitude Xo- within less 
than 1 %. Thus, irregularities in the amplitude are due 
not to a drift in time, but mainly to reversible effects of 
temperature variation. 

Putting aside the amplitude problem, our results are 
in excellent accord with the neutron scattering data of 
sect. II. This motivates us to reconsider dielectric loss 
data which have been analysed previously with A = 0.78. 

IV. OTHER SPECTROSCOPIES 
A. Dielectric Loss 

Dielectric spectroscopy on PC has been described re- 
cently j||2^,Q . The measurements extended over a wide 
range of temperatures and many decades in frequency, 
and the analysis has addressed different issues which are 
currently debated in the context of glass transition dy- 
namics. 



Within the GHz-THz range, the dielectric susceptibil- 
ity e"{uj) passes through a minimum as suggested by MC 
theory. Fits with Eq. (^), approximated as a sum of two 
power laws |^^, gave A = 0.78, and from the temperature 
dependence of frequency and amplitude (Q) Tc ~ 187 K 
was found. Though these values had strong support in 
the existing literature, they differ from our neutron and 
light scattering results. It is therefore interesting to ask 
whether the dielectric loss would also allow for A — 0.72. 

Since the available dielectric data are even noisier than 
the neutron susceptibilities, there is no unique way to de- 
termine the scaling range within which Eq. (|3|) applies; 
therefore, any value of A depends on the choice of the 
fit range. In the previous analysis |^^, fits were ap- 
plied between about 1 and 600 GHz. For temperatures 
around 200 K this range covered both sides of the mini- 
mum equally well. When A = 0.72 is imposed in an iter- 
ative master curve construction, the scaling range evades 
towards lower frequencies, as can be seen from the fits in 
Fig.|. 

For direct comparison. Fig. ^ shows master curves 
constructed with two different values of A. In the up- 
per curve, with A — 0.78, the measured susceptibilities 
have been rescaled with exactly the tOa and shown in 
Fig. 8 of Ref . . The rescaling works particularly well 
on the high-frequency side of the minimum and shows 
nicely the limit which has remained elusive in so many 
other experimental investigations. The lower curve, con- 
structed with A = 0.72, shows that the scaling property 
is completely lost above Wo-, whereas it has significantly 
improved in the low-frequency wing. Furthermore, the 
u!^^ limit of g\ describes the data nearly up to the a- 
relaxation maximum where they almost coincide, except 
for the lowest two temperatures. 

This unexpected observation motivates a new master 
curve construction which is quite common in the conven- 
tional analysis of a relaxation, but which has never before 
been applied to MC /3 relaxation: in Fig. |l^, the dielectric 
data are shown with original amplitude, rescaled only in 
frequency by the a maximum frequency Wmax- As usual 
for a relaxation (and in accord with the second scaling 
law of MC theory) the amplitude and the line shape of the 
a peak are temperature independent. But here, the scal- 
ing behaviour extends far down into the high-frequency 
wing — in fact, it extends as far as the data have been 
measured, except again the two lowest temperatures, 193 
and 203 K. 

However, a master curve which extends from the maxi- 
mum up to the minimum of a dynamic susceptibility can- 
not simultaneously obey the first and the second scaling 
law of MC: we expect 

Wmax/t^min K \a\^^^'' , (7) 

and, even more elementary, ernin oc |cr|"'^^^ whereas emax 
remains constant. 
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Since there is no doubt about the scahng of a relax- 
ation, we are bound to conclude that the dielectric data 
do not reach the asymptotic regime of fast /3 relaxation, 
except in a rather small temperature range that extends 
at best to about 210 K. For higher temperatures, the 
high-frequency wing of a relaxation, though technically 
describable by g\ , does not represent the first scaling law 
limit. 

This conclusion is independent of any fitting, and does 
in particular not depend upon an imposed value of A. 

B. Solvation Dynamics 

The solvation response of s-tetrazine has been mea- 
sured in PC from 1.5 to 100 ps A MC analysis in 
the time domain identified both a- and /^-scaling regions 
in the dynamics p6| . A unified analysis of both regions 
was consistent with MC theory for a range of A-Tc pairs. 
In the original publication, A = 0.78 was fixed, in ac- 
cord with the light scattering analysis of Du et al. . 
However, this value of A yielded a cross-over temperature 
Tc = 176 K substantially below the value obtained by Du 
et al. (Tc = 187 K). 

The solvation data have been reanalysed by the same 
method used in Ref. but using A = 0.72. The analy- 
sis of the ck-scaling region is unchanged. Figure ^ shows 
the new beta-scaling plot and the fit to g\{t/ta-), both 
of which are as good as in the previous analysis. The 
inset shows a temperature scaling plot of and r^. As 
expected, the scaling law deteriorates at temperatures 
far above Tc. Depending on the range of temperatures 
fit, acceptable values of Tc lie in the range 178-182 K. 
Thus, the solvation response is consistent with both the 
A and Tc obtained in Sects. II and III by neutron and 
light scattering. 

V. COMPARISON OF DYNAMIC OBSERVABLES 
A. Direct Comparison of Susceptibilities 

Neutron scattering, light scattering with or without 
proper bandpass, dielectric loss: each experiment, taken 
alone, seemed in full accord with the asymptotic predic- 
tions of MC theory. Taken together, the situation be- 
comes more complicated. 

In the fast (3 regime, any observable that couples to 
density fluctuations is expected to tend towards the same 
asymptotic limit (^). For neutron scattering, this pre- 
diction takes the form of a q,uj factorisation, and is con- 
firmed, in PC as in other materials, over a wide frequency 
range. However, it breaks down completely when depo- 
larised light scattering or dielectric loss are included. 

In Fig. dynamic susceptibilities from neutron scat- 
tering, light scattering [p9|, and dielectric loss are shown 



for direct comparison on an arbitrary intensity scale, but 
in absolute frequency units. The result is in flagrant con- 
tradiction to any factorisation property: there is no one 
temperature for which the measured e"(w) and fall 
together; in particular, their minimum frequencies differ 
systematically. 

On a first sight, this outcome is a bad surprise, and 
could make us doubt whether fits of individual data sets 
with the asymptotic laws and are meaningful at 
all. On closer examination, the discrepancies between 
the three data sets can be traced back to two major dif- 
ferences: to the individually different temperature ranges 
within which the /3 asymptote applies, and to a system- 
atic shift in the frequency scale lo^. 

B. Temperature Range of /3 Relaxation 

Within idealised MC theory, all dynamic observables 
converge towards the same scaling limit (||) and , char- 
acterised by just one lineshape parameter A and one fre- 
quency scale Lo^. This universality, however, is restricted 
to the lowest-order asymptote and does not imply a uni- 
versal radius of convergence: the next-to-leading-order 
corrections already depend on the microscopic coupling 
]l]| so that different observables may reach the asymp- 
tote at different temperatures and frequencies. 

Nevertheless, for neutron and light scattering, as well 
as for solvation response, (||) and (Q) hold over rather 
large temperature intervals. In neutron scattering, (^) 
holds best between 220 and 260 K; at 210 K, the suscep- 
tibility minimum approaches the instrumental resolution 
function, and the signal becomes very weak: these tech- 
nical limitations prevent us from following S{q^Lo) closer 
to Tc. For light scattering. Fig. || confirms (Q) with rea- 
sonable precision up to nearly 1.5 Tc. 

On the other hand, the a master curve in Fig. ^ 
suggested that the range of f3 scaling is rather small 
for dielectric loss. This is corroborated by Fig. |lj in 
which the uj^ from fits with fixed A = 0.72 are compiled. 
For the lowest temperatures 193-213 K, the minima of 
e"{iij) coincide with the light scattering data, whereas 
for higher temperatures the dielectric lu^ cross over to a 
much steeper temperature dependence characteristic for 
a relaxation (actually even steeper than the MC predic- 
tion \a\ ' ' j. 

C. Frequency Scales tUcr 

As Fig. ^ shows, the tUa- from neutron scattering fall 
about 35 % below the uj^ from light scattering. The good 
agreement of uj^ from neutron scattering with t^^^ from 
solvation dynamics excludes the possibility that the dis- 
crepancy between neutron and light scattering lies only 
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in multiple scattering or other technical shortcomings of 
neutron scattering. Ergo, at least one of the two scatter- 
ing techniques does not see the true /3 asymptote of MC 
theory. 

On the other hand, the observed scahng, the quality of 
the g\ fits, and the accord of the A and Tc call for a MC 
interpretation, and indicate a transient behaviour that is 
closely coupled to the universal (3 asymptote. 



VI. DISCUSSION 
A. Strength of a Relaxation 

Although there is no universal criterion for the validity 
of /3 scaling, an upper limit may be given: it is plausible 
that the asymptote (||) which is based on an expansion 
around the susceptibility minimum will no longer apply 
when the minimum ceases to exist. This happens when 
the height of the minimum, Xa-, becomes comparable to 
the maximum of x"{^) in the vibrational band. This con- 
dition also determines the highest temperature for which 
the proportionality (Q) can hold. 

On the other hand, the temperature range over 
which holds can be related to the strength of the 
a peak. To explain this relation, we refer back to the a 
master plot in Fig. 11, Besides the dielectric e" {ui / ujmayd > 
the figure also shows Xvni^ / ^max) from depolarised light 
scattering. For both data sets, the same frequency scale 
cjiiiax has been applied, as determined from the maximum 
position of the dielectric a peak. The intensity scale of 
the light scattering data is arbitrary and has been ad- 
justed by a global factor so that the maximum of the 
microscopic excitations has about the same height as for 
the dielectric data which are measured in absolute units. 

Compared to the high-frequency maximum, the a peak 
is about four times weaker in light scattering than in 
dielectric loss. In the high-frequency wing, at, say, 
uj ~ SOoJinax, the ratio is reduced to a factor of about 
2. This reduction is due to different positions of the a 
maxima, and to different slopes of the wings. 

At still higher frequencies, the susceptibilities cross 
over towards the minimum. For the lowest temperatures 
shown in the figure, the ratio of 2 between dielectric and 
light scattering data survives in this frequency range. 
Thus, for a low-temperature value of cr, the P ampli- 
tude XiT for light scattering is about two times smaller 
than for dielectric loss. 

On the other, from the condition stated above, and be- 
cause we have rescaled the susceptibilities to about the 
same values at high frequencies, the temperature range of 
f3 scaling ends for both techniques at about the same ab- 

1/2 

solute value of Xcr • Using the proportionality Xo- oc |(t| , 
the dielectric data reach this limit after a 4 times smaller 
temperature change than the light scattering data do. 



Thus, the strength of a relaxation can explain why the 
(3 scaling regime extends over only 30 K in dielectric loss, 
compared to 90 K or more in light scattering. 



B. Frequency Range of (3 Relaxation 

One can imagine many experimental imperfections 
that distort spectra measured in a scattering experiment. 
For instance, multiple scattering could overlay S{q, lo) 
with a convolution with itself. But any such distortions 
would affect the wings of the susceptibility much more 
than the region around the minimum. In particular, con- 
volutions lead to similar corrections as intrinsic next-to- 
leading-order terms. Therefore, it is not easy to explain 
the discrepancy observed in the ujcr- 

We note, however, that the gx fits to the neutron scat- 
tering data in Fig. |4| extend up to frequencies between 200 
and 400 GHz, whereas in light scattering, the fitted g\ al- 
ready deviate from the measured susceptibilities a little 
above 100 GHz. It is therefore quite possible that our 
neutron scattering analysis sees a preasymptotic tran- 
sient rather than the true gx. 

This would ressemble the situation in glycerol pO[ , 
where the x"i'-^) allowed for the construction of a master 
curve and followed some MC predictions, although the 
full f3 asymptote was not reached. Numeric solutions of 
a simple two-correlator model showed how such a sce- 
nario can arise from a MC ansatz ^l|. Similar fits, with 
one slave correlator for each observable, also work for 
propylene carbonate l42[ . 



C. Current Theoretical Developments 

Some of the questions raised by our experiments are 
also addressed by recent numerical solutions of MC equa- 
tions and molecular-dynamics simulations. 

A MC analysis of the hard-sphere liquid as a function 
of wavenumber has demonstrated that corrections to scal- 
ing are of differing importance for different observables. 
Analytic expansions have shown that corrections to the 
a process are of higher order than corrections to the /3 
asymptote: this may explain why a scaling holds over a 
wider range than the factorisation property of /3 relax- 
ation |l|]. 

In a molecular-dynamics simulation of a liquid made 
of rigid diatomic molecules , orientational correlation 
functions have been analysed for different angular mo- 
menta I — 0,1,2,... All correlators were found to fulfil 
the MC factorisation — with the pronounced exception 
of ^ = 1 for which the position of the susceptibility min- 
imum is shifted by a factor of 10 [loc. cit., Fig. 8]. The 
authors attribute this peculiar behaviour to 180° jumps, 
which is not necessarily the right explanation for similar 
behaviour in more complex fluids. A scenario with an 
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underlying type A transition for odd I has been proposed 
recently Q] . It has been noted before that the dielectric 
response couples only weakly to vibrational excitations 
, and the peculiar strength of the a peak in the I = 1 
correlator has also been found in an analytic extension 
of MC theory to non-spherical particles . 

A decomposition according to angular momenta might 
also explain the excellent accord between neutron scatter- 
ing and solvation response found in Fig. since both 
techniques see I — whereas light scattering might be 
dominated by ? = 2 [||| . 
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TABLE L Exponent parameter A and cross-over temper- 
ature Tc of propylene carbonate as determined by different 
experimental techniques. 



method 


par am. ^ 


X 


Tc (K) 


neutron scattering 








210 


Brillouin scattering 


h 






270 ±5 


Brillouin scattering 


h 




c 


c 


viscosity^ 


Ta 




0.70 


196 


same data 


Ta 




0.78'' 


188 ±3 


neutron scattering 


Ta 




0.70*" 


180/188 


dielectric loss 


Ta 




0.78'" 


187'' 


VH Ught scattering 


LOc 


± ' X(y± 


0.78 


187 ±5=* 


dielectric loss 


UJa 


, X<y 


0.78 


187 


solvation response" 


ilT ! 


Ta 


0.78'" 


176 


same data 


ta , 


Ta 


0.72'" 


180 ±2 


neutron scattering 


Ua 


, X<T 


0.72 


182 


VH light scattering 




, Xt 


0.72 


182 



reference 




this work 
this work 
this work 



''This column indicates whether the MC fit was based 
on the Debye- Waller factor (|^), the a time scale Ta 
(^, or the /3 relaxation parameters uJa- and Xo- (§)• The 
subscript a± indicates that T < data have also been 
used. 

''Used as input. 

'^Determination of MC parameters judged unreliable by 
the original authors. 

''The value 179 ± 2 K obtained in an alternative analysis 
using extended MC theory is not directly comparable to 
the Tc's from idealised theory. 

''Transient hole burning on s-tetrazine in propylene car- 
bonate. Original analysis. 
'Original data from Ref. |^. 



FIG. 1. Incoherent neutron scattering law S{q, cu) of propy- 
lene carbonate at T = 207 — 210 K and q = 1.2 A~', mea- 
sured on two time-of-flight spectrometers under different ex- 
perimental conditions: Mibemol (Ai = 8.5 A, circular Al con- 
tainer) and Neat (Ai = 5.5 A, linear arrangement of glass 
capillaries). Figure (b) repeats the data of (a) on a double 
logarithmic scale and compares them to the measured resolu- 
tion functions (lines). Above 20 GHz, the quasielastic spectra 
agree within about statistical error. 



FIG. 2. Decomposition of Xq{^) according to the factorisa- 
tion (|) for Mibemol (210 K) and Neat (207 K) data, (a) The 
amplitude hq has been determined by least-squares matching 
between 50 and 2500 GHz. Other methods yield almost iden- 
tical results. — (b) The spectral function x"{^) is obtained 
after dividing measured Xgi^) by hq. At low frequencies, the 
accord between both experiments is excellent; however, at 
high wavenumbers and especially at high frequencies, devia- 
tions are strong and systematic. Multiple scattering and im- 
perfect correction for energy-dependent detector efficiency are 
probably the main causes for distortions of the spectral line 
shape. 



FIG. 3. Rescaled susceptibilities Xq(uS)lhq. The resolu- 
tion-broadened elastic line has been cut off, except at 285 K 
where the g-dependent a relaxation can be fully resolved. The 
factorisation (n) is confirmed for both the mode-coupling /3 
regime and the vibrational spectra in the THz region. 



FIG. 4. Susceptibilities x"('^), obtained as g- independent 
averages x'q{'-^)/hq from Mibemol data, for temperatures vary- 
ing in 10 K steps from 210 to 260 K. Solid lines are fits with 
the mode coupling asymptote g\{Lo/LOa) (^) with A = 0.72 
fixed. The fit ranges, indicated by full symbols, have been 
determined in a self-consistent iteration from Fig. at low 
frequencies, they are limited by resolution or a relaxation, 
at high frequencies by the low-lying vibrational peak. — In 
the inset, the logarithms of the fit parameters Xt and ojo- are 
plotted against each other. The straight line shows the slope 
a = 0.318 expected from theory. 



FIG. 5. Susceptibility master curves x"('^)/Xct "^s. uj/ujc 
from the data set of Fig. |^. Fitting the same original data with 
different values of A yields different scales X" s-nd uJa which 
lead then to distinctly differing master curves. The temper- 
ature dependence (^ of X" (O) and ui^ (♦) is tested in the 
insets. — (a) Iterative rescaling yields a A of about 0.69 (but 
unstable and depending much upon idiosyncrasies of the fit- 
ting procedure). The temperature dependence of Xa a^nd uj^ 
is in conflict with (^). — (b) For A ~ 0.72, the rescaling is as 
good as for 0.69, and x<y and uja extrapolate from 220-260 K 
to a consistent Tc = 182 K. — (c) With A = 0.78 as suggested 
by part of the literature, the master curve is of poorer qual- 
ity, fits with g\ show that the imposed value of A is not self 
consistent, and X" and are again in confiict with (ku). 
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FIG. 6. Susceptibility xvn{^) from depolarised light scat- 
tering from propylene carbonate at 216 K, measured at the 
Fabry-Perot-Sandercock tandem interferometer at Garching, 
with Ai = 514.5 nm and with mirror spacings zo = 0.8,2.4 
and 7.2 mm. The matching of the three data sets is excel- 
lent and would extent even further without the VV Brillouin 
mode leaking around 10 GHz. The solid lines is a fit with the 
mode-coupling asymptote Xag\{uj/oJa) with A = 0.72. The 
dotted line corresponds to the mean literature value A = 0.78. 



FIG. 11. In this master plot, the frequency scale is taken 
from the position of the a relaxation peak in dielectric loss, 
and applied to both dielectric loss e"(<jj) (same data as 
and depolarised light scattering Xvh('^) (scanned from 
2l|]). While the e"((jj) are shown in absolute units, the 
amplitude of Xvh('^) is arbitrarily rescaled by a global fac- 
tor so that the microscopic peaks culminate at about equal 
height. — This scaling representation elucidates how a strong 
a peak reduces the temperature range for which a suscepti- 
bility minimum can be observed. 




FIG. 7. Depolarised light scattering susceptibilities, com- 
bined from two spectral ranges [zo = 0.8 and 2.4 mm), as 
function of temperature. Solid lines are mode-coupling fits 
with fixed A = 0.72. Open symbols indicate regions which 
had to be excluded from the fits; above about 80 GHz the 
cross-over to the non-universal vibrational spectrum; around 
8-15 GHz the Brillouin line, and at the highest temperature 
the beginning of a-relaxation. 

FIG. 8. Scaling parameters obtained from mode-coupling 
fits. The rectified plot (a) confirms that LJa evolves with 
When fitting only the fuU symbols, Tc = 182 ± IK is 
obtained in full accord with the neutron scattering result from 
Fig. 1^. When extending the temperature range, Tc tends to 
decrease by a few K. — The amplitudes, determined indepen- 
dently for two spectral ranges, are advantageously represented 
(b) in a logarithmic plot vs. T — Tc with Tc = 182 K as de- 
termined above. The solid line shows the slope 1/2 expected 
from theory. 

FIG. 9. Dielectric loss data as published in Refs. |p|, p3| , p4) . 
SoUd lines are MC fits (|) with fixed A = 0.72. This figure can 
be compared directly to Fig. 7 in Refs. where the same 
data are fitted with A = 0.78. 

FIG. 10. Dielectric loss data [^,^,^, rescaled to two ver- 
sions of a master function e' juja) j . In the upper curve, 
the same and are used as determined in Ref. for 
A = 0.78. For the lower curve, A = 0.72 was imposed which 
shifts the scaling range towards the low-frequency wing of the 
susceptibility minimum. 



FIG. 12. Solvation response of s-tetrazine in propylene car- 
bonate as measured by transient hole burning. Same data as 
in Fig. 6 of but now reduced with A = 0.72, /c = 0.56. 
The solid curve is the scaling function g\{t/tc). The devia- 
tions at long times are attributed to ct relaxation. — The inset 
(c/. Fig. 7 of j2^) shows the temperature scaling of ta {tp in 
Pg) and Ta (both in ps), using the exponents 2a — 0.636 
and 7 — 2.395 that correspond to A = 0.72. The fits shown 
extrapolate to a common Tc = 182 K in accord with neutron 
scattering measurements. (Good fits can be obtained over the 
range 178-182 K). 

FIG. 13. Direct comparison of susceptibilities from neu- 
tron scattering [Mibemol, this work], depolarised light scat- 
tering j2l], and dielectric loss [p|, p3|j24| ] . The absolute scale 
is chosen arbitrarily (but temperature-independent) to make 
the data approximately coincide in the low-temperature, 
low-frequency corner where they possibly reach a common 
asymptotic regime. 

FIG. 14. Frequencies cucr, determined from mode-coupling 
fits with the same A — 0.72, for neutron scattering, light 
scattering, dielectric loss, and solvation dynamics (from the 
time-domain analysis converted to o^o- = l/^o-). Temperatures 
are shown as T — Tc with Tc = 182 K. The solid line shows 
the slope l/2o = 1.572 expected for the given value of A; the 
dotted line, fitted to some dielectric data points, has a slope 
of 3.1 which is even steeper than l/2a + 1/26 — 2.4 expected 
for the a-relaxation limit. 
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